library("palmerpenguins")
library("ggplot2")
library("dplyr")
library("readxl")
bovine <- read_xlsx("data/bovine-tb/gb-tb-stats.xlsx") |>
mutate(date = as.Date(date), year = format(date, "%Y"),
doy = as.numeric(format(date, "%j"))) |>
rename(n_cases = n_not_otf)
penguin_labs <- labs(x = "Bill length (mm)", y = "Flipper length(mm)")